Electronic Properties of the Hubbard Model on a Frustrated Triangular Lattice 
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We study electronic properties of the Hubbard model on a triangular lattice using the cellular 
dynamical mean-field theory. The interplay of strong geometric frustration and electron correlations 
causes a Mott transition at the Hubbard interaction U/t = 10.5 and an unusual suppression of low 
energy spin excitations. Doping of a triangular Mott insulator leads to a quasiparticle peak (no 
pseudogap) at the Fermi surface and to an unexpected increase of low energy spin excitations, in 
stark contrast to the unfrustrated square lattice case. The present results give much insight into 
strongly frustrated electronic systems. A few predictions are made. 
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Geometric frustration with strong electronic correla- 
tions is one of the main issues in modern condensed mat- 
ter physics. The simplest example is the two dimensional 
Heisenberg model (large U limit of the half-filled Hub- 
bard model) on a triangular lattice in which all three 
spins cannot be antiparallel at the same time. The frus- 
trated triangular lattice geometry was argued by Ander- 
son [l[ to provide an ideal background for the long sought 
resonating valence bond (RVB) state. The interplay of 
strong geometric frustration and electronic correlations is 
expected to lead to some exotic phases. Recent discovery 
of superconductivity in the triangular lattice compound 
Na^CoO^ • yH^O |2| and a possible quantum spin liquid 
state in the anisotropic triangular lattice organic mate- 
rial K -(ET) 2 Cu 2 (CN) 3 1 has further stimulated interest 
in strongly frustrated electronic systems. 

The electronic properties of the Hubbard and t — J 
models on a triangular lattice have been studied using 
various analytical and numerical techniques. These in- 
clude a high-temperature expansion [4], a slave-boson 
mean- field [a, la 0, an RG method an exact diag- 
onalization (ED) technique fiol ITll ] , a quantum Monte 
Carlo (QMC) simulation [Tj|], a dynamical mean-field 
theory (DMFT) HE03, and a variational Monte Carlo 
approach [l5j . However, some of the central issues on 
a triangular lattice have not been clearly addressed and 
still remain an open question: what is the nature of the 
ground state of a doped triangular Mott insulator? is 
it a non-Fermi liquid (manifested as the presence of a 
pseudogap at small doping like on an unfrustrated square 
lattice) or a correlated Fermi liquid? what are qualita- 
tive differences between the triangular and square lat- 
tice systems? In this paper, utilizing recent theoretical 
progress in quantum cluster methods, we address these 
issues and provide much insight into electronic properties 
in strongly frustrated electronic systems. 

The two-dimensional Hubbard model on a triangular 
lattice is described by 
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FIG. 1: (a) Three-site quantum cluster embedded in an effec- 
tive medium (called bath) and (b) first Brillouin zone of the 
triangular lattice. 



electrons of spin a, — (^ i<T c%a is the density of a spin 
electrons, U is the on-site repulsive interaction, and \± is 
the chemical potential controlling the electron density. 

To describe strong electron correlations and geometric 
frustration accurately, we use cellular dynamical mean- 
field theory (CDMFT) [l6[, a quantum cluster approach 
that allows one to extend DMFT 17j to incorporate 
short-range spatial correlations explicitly. CDMFT has 
been benchmarked and is accurate even in one dimen- 
sion nam. The infinite lattice is tiled with identical 
clusters of size N c , and the degrees of freedom in the 
cluster are treated exactly while the remaining ones are 
replaced by a bath of non-interacting electrons that is de- 
termined self-consistently (Fig. HJa)). To solve the quan- 
tum cluster embedded in an effective medium, we con- 
sider a cluster-bath Hamiltonian of the form [la HI, H3| 
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where cL 



(cicr) are creation (annihilation) operators for 



Here the indices fi, v — 1, • • • , N c label sites within the 
cluster, and and a^ a annihilate electrons on the clus- 
ter and the bath, respectively, are the hopping matrix 
elements within the cluster (the chemical potential fj, is 
absorbed here), £^ CT are the bath energies and V^„ CT are 
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FIG. 2: Density of states N(u>) on a triangular lattice at half- 
filling (n — 1) with increasing U/t. 



the bath-cluster hybridization matrices. The exact diag- 
onalization method [2l[ is used to solve the cluster-bath 
Hamiltonian Eq. [5] at zero temperature, which has the 
advantage of computing dynamical quantities directly in 
real frequency and of treating the large U regime without 
difficulty. In the present study we used N c = 3 sites for 
the cluster and Nb — 9 sites for the bath with m = 1, 2, 3, 
a = 1,2,3. Although the present study focuses on a small 
cluster with additional 9 bath sites, we expect our results 
to be robust with respect to an increase in the cluster 
size. This was verified by our recent low (but finite) 
temperature CDMFT+QMC calculations []J| where at 
intermediate to strong coupling a 4-site cluster accounts 
for more than 95% of the correlation effect of the infinite 
size cluster in the single-particle spectrum. Quantita- 
tively similar results were found for different bath sizes 
(Nb = 6). In this paper we focus on the normal state 
properties at and near half-filling. 

We first study the Mott transition of the half-filled 
Hubbard model on a triangular lattice. Among several 
physical quantities one can gauge to investigate the Mott 
transition, we present the density of states N(to) as a 
function of U/t shown in Fig. [2j For U/t < 10.5 a small 
spectral weight remains near the Fermi energy reminis- 
cent of the Kondo resonance observed in DMFT p7| . 
When U/t approaches 10.5, the spectral weight of the 
central peak is completely pushed away from u> = lead- 
ing to an interaction-driven Mott insulator. This critical 
value of U/t is somewhat smaller than 12 obtained in 
the ED study of an isolated 12-site cluster [9| and in the 
DMFT analysis [II] . Beyond U/t = 10.5 it becomes in- 
creasingly more difficult to find a convergent solution, 
since the ground state starts to become degenerate. The 
way how the Mott transition occurs on a frustrated tri- 
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FIG. 3: (color online) Imaginary part of the local dynamical 
(a) charge and (b) spin correlation functions at half-filling 
with increasing U/t. The dashed, dotted and solid curves 
correspond to U /t=9, 10 and 10.5, respectively. 



angular lattice is dramatically different from the unfrus- 
trated square lattice case. In the latter [2Cj], with in- 
creasing U/t the spectral weight always becomes a local 
minimum at the Fermi energy until it completely van- 
ishes near U/t = 6. 

Next we examine the local dynamical charge and spin 
correlation functions near the Mott transition shown in 
Fig. [3J With increasing U/t the low energy charge ex- 
citations are gradually depleted and transfered to high 
energy near lu — U/t. On the other hand, the low energy 
spectrum in the spin correlation function undergoes a 
dramatic change near the Mott transition. As the Mott 
transition is approached, the primary low energy exci- 
tation near w/t — 0.2 moves rapidly to the secondary 
peak near w/t = 0.5 — 0.6 and eventually disappears. 
Namely, geometric frustration suppresses the low energy 
spin excitations, which is most strongly manifested in 
the insulating state (U/t = 10.5). We found similar re- 
sults for the nearest neighbor charge and spin correlation 
functions. This feature is in stark contrast to the square 
lattice case. In the latter, for the same range of U/t 
(9 — 10.5) the low energy spin spectra are almost identi- 
cal with a sharp primary peak near w/t = 0.25 — 0.3 [22| 
and a secondary peak near w/t — 0.6 — 0.7 similar to the 
case of U/t = 9 on a triangular lattice. Near the range of 
U/t (4 — 6) where a continuous crossover of a metal to an 
insulator occurs on a square lattice, with increasing U/t 
the low energy spin spectrum near w/t = 0.25 becomes 
stronger, which is opposite to the behavior found on a 
triangular lattice. 

Figure 2] shows the spectral function A(k, lS) for both 
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FIG. 4: Spectral function A(k,u>) for (a) n=0.9, (b) n=l, and 
(c) n=l.l along the path r-M-K-F for U/t = 10.5. F, M, and 
K are defined in Fig. [lib). 



10% hole- and electron-dopings. Electron-doping corre- 
sponds to the negative sign of the hopping integral (t < 0) 
in the t — J model. At half-filling (Fig. g[b)) A(k,w) 
shows several features, which are similar to those on a 
square lattice [20(. Although they are weaker than on a 
square lattice, the low energy inner bands are clearly seen 
near the Fermi e nerg y, which are caused by short-range 
spin correlations [20(. Upon hole- or electron-doping the 
Fermi level jumps to the nearest low energy band just 
like on a square lattice. This indicates the importance of 
short-range correlations even on a highly frustrated trian- 
gular lattice. The most dramatic difference, however, is 
that near half-filling there is no evidence of a (strong cou- 
pling) pseudogap unlike on a square lattice. It is more 
evident from the fact that A(k,oj) becomes sharper as 
k approaches the Fermi surface. This is true for other 
fillings (5 — 30% hole- and electron-doping). Thus the 
present result clearly shows that a doped triangular Mott 
insulator is a correlated Fermi liquid unlike on a square 
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FIG. 5: (color online) Imaginary part of the local dynamical 
(a) charge and (b) spin correlation functions for U/t = 10.5 
with doping. The inset in (b) is the imaginary part of the 
local dynamical spin correlation function for the unfrustrated 
square lattice case. The dashed, solid, and dotted curves 
correspond to n=Q.9, 1 and 1.1, respectively. 



lattice [20| ■ The absence of the pseudogap is ascribed to 
too much frustration of short-range spin correlations on a 
triangular lattice as shown in Fig.[BJ Here we predict that 
near half-filling A(k, u>) has a quasiparticle peak at the 
Fermi surface, which should be tested by angle resolved 
photoemission spectroscopy (ARPES) experiments on a 
triangular lattice compound such as Na^OoC^ at small 
doping. 

Next let us examine the local dynamical charge and 
spin correlation functions upon doping shown in Fig. [5J 
As expected, with doping the high energy charge exci- 
tations at half-filling move to the low energy due to the 
metallic nature of electrons. The most surprising result 
comes from the spin correlation function. The primary 
low energy peak near w/t = 0.2, which has vanished in 
the insulating state due to strong geometric frustration 
shown in Fig. [3j reappears upon doping. This is com- 
pletely opposite to the square lattice case (the inset in 
Fig. EJb)) where the primary low energy peak weakens 
rapidly with doping. On a triangular lattice geometric 
frustration is released by doping. Here we also predict 
that inelastic neutron scattering experiments should ob- 
serve the increase of low energy spin excitations with 
doping if the undoped system is a Mott insulator (see 
caveat [Hj|). 

Finally we study the static local and nearest neigh- 
bor spin correlations with doping for both triangular and 
square lattices shown in Fig. [5] For the local quantity 
(Sf ■ Sf) they are similar in both cases, although the 
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FIG. 6: (color online) Static (a) local and (b) nearest neighbor 
spin correlations for U/t — 10.5 with doping. The dashed 
curves are the corresponding results on a square lattice. 

particle-hole symmetry is broken on a triangular lattice. 
The most dramatic difference is seen in the nearest neigh- 
bor spin correlation (Sf ■ First of all the nearest 
neighbor spin correlation is much reduced compared with 
a square lattice because of strong geometric frustration. 
While it remains essentially antiferromagnetic (negative 
sign) on a square lattice, it becomes rapidly ferromag- 
netic on the electron-doped side of the triangular lattice. 
The strong asymmetric behavior on a triangular lattice 
that the spin correlation becomes ferromagnetic only on 
the electron-doping side is in agreement with the high 
temperature expansion study of the t — J model by Ko- 
retsune et al. [J] and the analysis of an isolated three-site 
cluster by Merino et al. [l3|. It is also noteworthy that 
the ferromagnetic correlation is maximum near n = 1.35 
(x = 0.35), where superconductivity has been observed 
in Naa;Co02 • yEbO [2]. Although the small cluster used 
in this work (N c = 3) does not allow us to study which 
gap symmetry (p or f) is more stable and how robust 
this result would be with respect to several parameters 



such as doping level, our present result offers a natural 
route to spin triplet superconductivity in this compound 
within the one-band Hubbard model (without assuming 
the existence of unobserved hole pockets near k = K). 

Sodium cobaltatc Na x Co02 has shown a rich phase di- 
agram [24| . These include a superconducting phase near 
x = 0.3 — 0.35 upon hydration, a charge-ordered insu- 
lating phase at x = 0.5 and a spin-density-wave metallic 
phase above x = 0.75. Since the charge-ordered insulat- 
ing phase in cobaltates is observed at x = 1/2 doping 
which is not a natural commensurate value (x = 1/3 
or 2/3) of the triangular lattice, the one-band Hubbard 
model alone would not be enough to describe this phase. 
It requires longer-range Coulomb interaction Vij, which 
we did not consider in this work. 

In conclusion, we have studied electronic properties of 
the Hubbard model on a triangular lattice using the cellu- 
lar dynamical mean-field theory. The interplay of strong 
geometric frustration and electron correlations is respon- 
sible for many unconventional features, which are in stark 
contrast to the unfrustrated square lattice case. In par- 
ticular the appearance of a Fermi liquid upon doping is 
ascribed to too much frustration of short-range spin cor- 
relations on a triangular lattice. We predict that upon 
doping a triangular Mott insulator ARPES should ob- 
serve a quasiparticle peak (no pseudogap) at the Fermi 
surface and that inelastic neutron scattering should find 
an increase of low energy spin excitations. 
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